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FOREWORD 


This  reprint  report  comprises  two  parts.  The  first  part  (beginning  on  pg.  1)  reproduces  a 
paper  presented  at  the  7th  European  Conference  on  Underwater  Acoustics  on  5-8  July  2004.  The 
second  part  (beginning  on  pg.  8)  reproduces  a  related  presentation  made  at  the  same  conference. 
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A  comparison  of  the  passive  ranging  performance  between  two  hypothetical  sonar  systems  is 
presented  using  theoretical  methods.  Each  system  is  composed  of  eight  hydrophones  that  are 
symmetrically  located  in  a  straight  line  spanning  the  same  total  array  length.  The  systems  differ 
in  the  way  the  hydrophones  have  been  grouped,  which  significantly  affects  the  associated  range- 
estimator  variance.  The  performance  of  each  system  has  been  analyzed  using  triangulation- 
based  ranging  as  well  as  optimal  processing  that  would  asymptotically  achieve  the  Cramer-Rao 
lower  bound  (CRLB).  While  this  work  compares  two  specific  eight-sensor  array  configurations, 
the  results  have  implications  regarding  the  design  and  performance  of  any  multi-aperture  sonar 
system. 


1.  INTRODUCTION 

The  localization  performance  of  a  passive  sonar  system  varies  with  the  placement  of  the  sensors 
comprising  the  sonar  array.  In  many  applications,  the  accurate  estimation  of  range  from  an  array 
to  a  radiating  source  relies  on  the  system’s  ability  to  measure  wave-front  curvature.  As  a  result, 
configurations  designed  optimally  for  ranging  have  a  portion  of  the  sensors  placed  in  the  middle  of 
the  array  and  others  placed  near  the  ends.  Carter  [1]  showed,  for  a  symmetric  line  array,  that  the 
variance  of  the  maximum  likelihood  (ML)  range  estimate  is  minimized  when  the  sensors  are 
configured  in  groups  of  one  quarter,  one  half,  and  one  quarter.  Other  research  and  studies, 
consistent  with  this  approach,  supported  the  decision  to  install  triple-aperture  sonars  on  modem 
naval  platforms.  While  performance  has  been  proven  both  theoretically  and  at-sea,  a  three-site 
installation  is  not  always  economically  feasible.  Whether  the  sensors  are  mounted  on  the  hull  of  a 
submarine,  for  example,  or  moored  on  the  ocean  floor,  cost  savings  could  be  realized  by 
eliminating  the  middle  sensor  grouping.  The  objective  of  this  work  is  to  quantify  the  theoretical 
loss  in  ranging  performance  due  to  the  installation  of  two  groups  of  sensors  instead  of  three. 
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Two  hypothetical  sonar  systems  based  on  symmetric  line  arrays  have  been  used  to  perform 
this  analysis.  Each  system  is  composed  of  eight  sensors  that  span  the  same  overall  array  length 
and  differ  only  in  the  way  the  sensors  are  grouped.  The  performance  of  each  system  has  been 
analyzed  using  triangulation-based  ranging  as  well  as  optimum  processing  that  would 
asymptotically  achieve  the  Cramer-Rao  lower  bound  (CRLB)  [2].  The  choice  of  eight  sensors 
allows  for  a  comparison  of  systems  with  the  same  number  of  hydrophones  and  opportunity  for 
either  system  to  achieve  the  same  maximize  output  signal-to-noise  ratio  (SNR).  However,  there 
is  nothing  special  about  the  number  eight,  and  this  analysis  generally  applies  to  systems  with  an 
arbitrary  number  of  sensors. 

Fig.  1(a)  illustrates  the  first  sonar  system,  designated  System  I,  which  has  its  sensors  grouped 
in  two,  four-element  sub-arrays,  each  spanning  length,  Ls.  It  has  been  assumed  that  this  sub¬ 
array  length  is  much  smaller  than  the  total  array  length,  or  L»Ls.  The  range  and  bearing  to  the 
source  referenced  to  the  central  origin  is  denoted  by  R  and  B,  respectively.  Using  triangulation, 


Fig.l:  (a)  System  I,  an  eight-sensor,  dual-aperture  configuration,  (b)  System  II,  a 
triple-aperture  configuration  formed  by  grouping  four  sensors  in  the  middle  of  the  array. 

range  R  can  be  expressed  in  terms  of  the  System-I  sub-array  bearings,  Bn  and  Bn,  which  are 
approximately  separated  by  the  total  array  length,  L,  given  that  L»Ls.  Fig.  1(b)  illustrates  the 
second  sonar  system,  designated  System  II,  which  has  its  sensors  configured  optimally  in  groups 
of  one  quarter,  one  half,  and  one  quarter.  The  middle  grouping  of  four  sensors  spans  the  same 
length,  Ls,  as  each  of  the  sub-arrays  in  System  I.  However,  the  System-II  sub-arrays  formed 
using  two  middle  sensors  and  two  end  sensors  span  half  the  total  array  length  and  have  phase 
centers  located  at  x  =  ±L/ 4.  Thus,  range  can  be  expressed  in  terms  of  the  bearings,  Bn\  and  Bin, 
which  are  approximately  separated  by  LI 2. 

The  analysis  that  follows  assumes  what  Van  Trees  [2]  refers  to  as  SPLOT,  i.e.,  stationary 
processes,  long  observation  times.  Furthermore,  signals  and  noises  are  assumed  Gaussian  and 
mutually  uncorrelated  with  the  same  arbitrary  power  spectra  at  each  sensor.  Thus,  a  direct-path. 
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isotropic-noise  model  is  implied,  where  the  propagated  acoustic  energy  to  each  sensor  is  the 
same.  There  is  also  an  implicit  assumption  that  the  SNR  is  high  enough  for  the  system  to  be 
tracking  the  local  peak  in  the  correlation  domain  and,  thereby,  achieve  performance  that  depends 
on  the  bandwidth  and  the  square  of  the  center  frequency.  This  is  classical  CRLB  theory  that  is 
difficult  to  achieve  below  a  threshold  SNR,  where  a  tighter  bound,  the  Ziv-Zakai  lower  bound 
(ZZLB)  applies.  Finally,  there  is  an  assumption  that  for  the  integration  time,  T,  the  source  and 
receiver  motion  is  negligible. 


2.  OPTIMUM  PERFORMANCE 


A  system’s  performance  is  commonly  measured  by  the  mean-squared  error  (MSE)  of  the 
parameter,  or  parameters,  it  is  used  to  estimate.  Optimum  performance  in  this  sense  is  achieved 
when  the  MSE  of  the  parameter  estimate  equals  the  Cramer-Rao  lower  bound  (CRLB)  [2].  If  an 
estimator  that  achieves  the  CRLB  exists,  it  is  the  maximum  likelihood  (ML)  estimator.  The  ML 
estimator  is  asymptotically  efficient ;  that  is,  its  mean-squared  error  approaches  the  CRLB  as 
more  data  are  incorporated  into  its  solution.  Hence,  for  SPLOT  conditions,  ML  estimates  of 
either  range  or  bearing  may  be  considered  optimum  in  the  MSE  sense.  The  statistics  of  the  ML 
estimator  were  derived  by  Bangs  and  Schultheiss  [3].  Under  their  assumptions  of  Gaussian 
random  processes  and  spatially  incoherent  noise,  it  was  shown  that  this  estimator  is  unbiased, 
and  thus  MSE  is  equivalent  to  the  variance  of  the  parameter.  Letting  £  denote  the  parameter  of 
interest  (range  or  bearing),  the  variance  of  the  ML  estimate  for  an  M- sensor  array  with  generic 
configuration  “A”,  is  given  by  the  following  [3]: 


JT_ 
2  n 


P '>) 

\  +  M-P(co) 


-1 


(1) 


In  Eq.  (1 ),  P(a> )  denotes  the  frequency  dependent  SNR,  S(a>)/N(oj),  where  S(co)  and  N(a>) 
represent  the  known  power  spectra  of  the  signal  and  noise,  respectively,  and  co  is  radian 
frequency.  M  is  equal  to  the  number  of  sensors  or  hydrophones  in  the  array,  Q  is  the  processing 
bandwidth  in  radians  per  second,  and  T  is  the  observation  time  in  seconds.  K(P  depends  on  the 
parameter  being  estimated,  as  well  as  the  placement  of  the  array’s  sensors.  Let£  =  R ,  where  R 
denotes  the  range  from  the  origin  to  the  source,  and  let  B  denote  the  bearing  angle  relative  to  the 
array’s  axis,  as  illustrated  for  both  systems  in  Fig.l.  Eq.  (1)  represents  the  range  variance  when 
KpR)  is  given  by  [3] 
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where  xm  is  the  coordinate  of  sensor  m  along  the  axis  of  the  array  in  meters,  and  c  is  the  speed  of 
sound  in  water  in  meters  per  second.  Eq.  (2)  is  valid  when  the  sensors  are  placed  symmetrically 
about  the  array’s  origin,  which  is  consistent  with  the  analysis  presented  here.  Analyses  of 
asymmetric  arrays  with  arbitrary  origin  placement  require  an  additional  degradation  factor. 
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which  would  account  for  the  coupling  of  range  and  bearing  uncertainty  [3].  Eqs.  (1)  and  (2)  hold 
only  for  the  estimation  of  a  single  parameter  (range),  given  that  all  other  parameters  are  known. 
Note  that  the  contribution  of  Eq.  (2)  to  the  range  variance  is  purely  geometrical.  Thus,  Eq.  (1)  is 
the  product  of  two  distinct  components:  a  signal  dependent  factor  and  a  geometrically  dependent 
factor.  The  former  includes  quantities  such  as  frequency-dependent  SNRs  at  the  sensors, 
bandwidth,  and  observation  time.  Thus,  if  the  same  optimum  processing  is  applied  to  two 
Af-sensor  array  configurations  under  the  same  signal  and  noise  conditions,  then  the  difference  in 
ranging  performance  depends  only  on  the  placement  of  the  array’s  sensors. 

Eq.  (2)  can  be  computed  for  each  sonar  system  and  substituted  into  Eq.  (1)  to  yield  the  range 
variance  in  terms  of  sub-array  length,  4s,  and  total  array  length,  L.  For  System  I  (A=l)  and 
System  II  ( A=I1 ),  respectively, 


4L\l)smaB  \2n  J°  l+8/>(®)J 
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2 c2R4  f  T  rn  2  p2{(o ) 
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(3) 

(4) 


Hence,  the  range  variance,  or  MSE,  for  System  II  is  independent  of  4s  under  the  assumption 
of  L»Ls  and  is  inversely  proportional  to  the  fourth  power  of  the  total  array  length,  L.  This  is 
consistent  with  Carter’s  result  in  [1]  for  the  optimal  ranging  array  with  quarter,  half,  quarter 
sensor  groupings.  In  contrast,  the  range  variance  for  the  dual-aperture  configuration  of  System  I 
is  inversely  proportional  to  the  squared  product  of  L  and  L$.  The  relative  performance  is 
expressed  by  taking  the  square  root  of  the  ratio  of  Eqs.  (3)  and  (4): 
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Eq.  (5)  shows  that  the  range  error  of  System  I  is  proportional  to  the  range  error  of  System  II 
and  grows  linearly  with  the  ratio  of  L  to  Ls.  When  this  ratio  is  equal  to  thirty,  for  example,  the 
range  error  for  the  dual-aperture  system,  System  I,  would  be  ten  times  the  error  expected  for  the 
triple-aperture  system.  System  II. 


3.  TRIANGULATION 

The  result  summarized  by  Eq.  (5)  states  that  System  II  theoretically  outperforms  System  I  if 
optimum  processing  were  attainable.  A  more  practical  comparison  would  be  one  that  accounts 
for  a  potential  loss  in  performance  due  to  sub-optimal  processing.  Such  a  comparison  can  be 
made  for  triangulation-based  ranging,  which  relies  on  independent  relative  bearing  estimates 
from  two  sub-arrays  and  a  simple  geometrical  model  to  triangulate  the  range  from  the  array  to 
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the  acoustic  source.  The  range  error  variance  due  to  triangulation  can  be  expressed 
approximately  by  [4]: 

4-  [  ^2(^i)+ ^AijX  (6) 

le 

where  cr2  (5^1)  and  g2[ba2J  are  the  variances  of  the  two  bearings  used  to  estimate  range,  R  is 
the  actual  range  to  the  source,  and  Le  is  the  effective  phase-center  separation  of  the  sub-arrays 
that  independently  produce  the  bearing  estimates.  System  I  generates  bearing 
estimates  Bn  and  Bl2  using  sub-arrays  that  are  effectively  separated  by  LE  =  L  sin  B ,  as  shown  by 
Fig.  1(a).  System  II  generates  bearing  estimates  Bin  and  Bin  using  sub-arrays  that  are 
effectively  separated  by  LE  -  L/2sin  B,  as  shown  by  Fig.  1(b).  Note  that  in  Eq.  (6),  and  in  what 
follows,  the  under-bar  notation,  A ,  identifies  that  two  independent  sub-array  processes  are 
performed  using  half  the  sensors  in  generic  array  configuration  “T”  to  produce  a  final  range 
estimate  via  triangulation. 

Triangulation  performance  can  be  related  to  the  underlying  array  properties  and  SNR 
conditions  using  Eq.  (1 )  to  obtain  the  ML  bearing  variance  for  the  sub-arrays  in  each  system. 

For  £  =  B ,  the  geometrically  dependent  factor  in  Eq.  (1)  for  an  M/2 -element  sub-array  is  defined 
as  follows  [3]: 


(7) 


Unlike  its  range  counterpart,  bearing  error  variance  is  virtually  independent  of  the  defined 
origin  since  Eq.  (7)  depends  on  the  differences  in  sensor  coordinates  instead  of  the  differences  of 
the  squared  coordinates  as  given  by  Eq.  (2).  Therefore,  there  is  no  need  to  redefine  the  origin 
placement  for  each  sub-array  analysis.  Using  the  appropriate  sensor  separations,  Eq.  (7)  can  be 
used  to  compute  K\B)  and  ,  the  geometrically  dependent  factors  of  the  sub-array  bearing 

variances  for  System  I  and  System  II,  respectively.  Due  to  the  symmetry  of  the  arrays, 

applies  to  both  sub-arrays  in  System  I,  and  KjP applies  to  both  sub-arrays  in  System  II. 

This  symmetry,  in  conjunction  with  the  assumptions  made  regarding  signal  and  noise  at  each 

sensor,  yields  equal  bearing  variances  for  the  sub-arrays  in  System  I, 

g2 (£/,  )=  g2 {Bn )=  g2 (£, ),  and  in  System  II,  g2 (£,/,)=  g2 [bi]2)  =  g2 [bi± ). 

By  substituting  the  appropriate  bearing  variances  and  effective  sub-array  separations,  as 
defined  above,  into  Eq.  (6),  the  triangulation  range  variance  for  each  system  is  obtained: 
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Comparing  Eqs.  (8)  and  (9)  with  Eqs.  (3)  and  (4)  shows  that  triangulation  with  ML  bearing 
estimates  theoretically  achieves  the  CRLB  when  the  output  SNR  is  high,  or  MP(co)»  1 .  When 
the  output  SNR  is  low,  the  efficiency  of  triangulation  degrades  by  seventy  percent  due  to  the 
independent  processing  of  each  M/2-element  sub-array  (the  triangulation  range  variance  is  twice 
the  CRLB).  However,  as  long  as  both  systems  operate  under  the  same  input  SNR  conditions, 
their  relative  triangulation  performance  is  equivalent  to  the  relative  performance  seen  for 
optimum  processing,  i.e.  the  ratio  of  Eqs.  (8)  and  (9)  is  equal  to  Eq.  (5).  This  is  an  intuitively 
satisfying  result:  System  II  has  an  inherent  ability  to  outperform  System  I  regardless  of  the  range 
estimation  methods  evaluated  here.  While  this  comparison  is  based  on  optimum  bearing 
estimation,  the  efficiency  of  triangulation  generally  scales  with  the  efficiency  of  the  underlying 
bearing  estimator.  It  can  be  shown  by  replacing  the  ML  bearing  estimator  with  a  split-beam 
tracker  [3],  for  example,  that  the  range  variances  given  by  Eqs.  (8)  and  (9)  both  increase  by  a 
factor  of  two.  The  resulting  ratio  comparing  System  I  to  System  II  remains  unchanged  and  is 
equal  to  Eq.  (5). 


4.  DISCUSSION  AND  CONCLUSIONS 


Under  the  assumptions  made  and  for  the  examples  cited,  it  has  been  shown  that  the  triple¬ 
aperture  system,  System  II,  shown  in  Fig.  1(b),  significantly  outperforms  the  dual-aperture 
system,  System  I,  shown  in  Fig.  1(a).  In  particular,  the  noise-induced  mean  squared  range  error 
of  System  I  can  be  reduced  by  a  factor  proportional  to  the  ratio  of  total  array  length  to  sub¬ 
aperture  length  by  simply  repositioning  half  the  sensors  in  the  middle  of  the  array,  as  in  System 
II.  While  this  comparison  was  performed  using  specific  eight-element  arrays,  the  results 
reported  here  apply  generally  to  any  two  M-element  symmetric  line  arrays  where  one  in 
configured  with  M/2,  M/2  sensor  groupings,  and  the  other  is  configured  with  M/4,  M/2,  M/4 
groupings.  In  fact,  by  extending  the  triangulation-based  analysis  in  Section  3,  the  ranging 
performance  comparison  can  be  generalized  to  yield  the  following: 


for  M  »  1 . 


(10) 


Eq.  (10)  is  very  similar  to  Eq.  (5);  both  are  linearly  related  to  the  ratio  of  L  to  Ls.  This  more 
general  result  shows  that  increasing  the  total  number  of  sensors  does  not  affect  relative 
performance.  Moreover,  for  dual-aperture  systems  with  sub-arrays  separated  by  tens  or  hundreds 
of  times  their  sub-aperture  lengths,  the  installation  of  a  third  sensor  site  could  reduce  the  range 
estimation  error  by  one  or  two  orders  of  magnitude.  Hence,  the  cost  saved  by  installing  two 
sensor  sites  instead  of  three  must  be  carefully  weighed  against  the  significantly  increased  ranging 
performance  of  a  three-site  system  composed  of  the  same  number  of  individual  hydrophone 
elements. 
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Problem  Statement  and  Assumptions 


>  Compare  the  underwater  acoustic  ranging  performance  of 
two  systems:  System  I  uses  dual  apertures;  System  II  uses 
triple  apertures.  (Keep  key  system,  processing,  and 
underwater  acoustic  parameters  constant.) 


>  Assumptions: 

•  Stationary  Processes,  Long  Observation  Times  (SPLOT) 

•  Signals  and  noises  are  Gaussian  and  mutually 
uncorrelated  with  the  same  arbitrary  power  spectra  at 
each  sensor 

•  Isotropic-noise,  direct-path  propagation,  propagated 
acoustic  energy  to  each  sensor  is  the  same  (SNR  at  the 
sensor  is  key  and  propagation  decouples  for  direct  path.) 
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Two  Sonar  Systems  Considered 


System  I  — q-q-q-q - * - o-o-o-o — 

Dual  Aperture  0  ^ 

s  s 

System  II  — <jx> - a-o+o-o - o-<j j> — ^ 

Triple  Aperture  ^ 

•  Two  eight-sensor,  symmetric  line  arrays  with  total  length  I 

•  System  I  is  a  dual-aperture  array  with  sub-aperture  length  Ls 

•  System  II  is  a  triple-aperture  array  which  is  formed  by 
grouping  four  sensors  in  the  middle  of  the  array  and  has  sub¬ 
aperture  lengths  Ls  and  approximately  LJ 2 
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The  Passive  Ranging  Geometry 


Source 


R  =  range  to  source  (meters) 

B  =  relative  bearing  to  source  (radians)  / 
Xj  =  coordinate  of  sensor  /  (meters) 

M-  number  of  sensors 


/  Range:  R 


A  Bearing:  B 


X,  x2  x3  ... 
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Qjjggj  Minimum  MSE  Ranging  Performance*  Saofe 


°jt  (*  )  =  KiR)  M 


Geometrical  Factor 


2  M  2  -  P 2  (a;) 
1  +  M-P(o) ) 


Signal-Dependent  Factor  increasing 


•  Depends  on  relative 
placement  of  the  array’s 
sensors  and  acoustic 
source 

•K  is  proportional  to  1/M2 

Decreasing  K  is  good 

*See  references,  including  paper  by 
Bangs  and  Scbullheiss  [1],  1973 


Where  T  and  M  is  good 

T  =  observation  time  (seconds) 

M  =  number  of  sensors 

to  =  frequency  (radians/sec) 

Q  =  bandwidth  (radians/sec) 

P(a>)  =  S(ro)/N(ro),  SNR  vs.  co 
S(a>)  =  signal  power  at  each  sensor 
N(cd)  =  noise  power  at  each  sensor 


July  5-8,  2004  Seventh  European  Conference  on  Underwater  Acoustics 

For  more  information  contac:  DiBiaseJH@NDt.NUWC.Navv.Mil  CarterGC@Npt.NUWC.Naw  Mil 


Components  of  the  Geometrical  Factor 


^=4c2^-  XI  (2-4) 


M  M 


(Note  K  is 
proportional 
to  1/M2) 


Sr-  sin  B  I  M  j=l 

\ _ t _ )  v  J 


-1 


Depends  on  sound 
speed,  c  (meters/sec) 


Depends  on  array  configuration,  i.e., 
placement  of  the  sensors 


Depends  on 
source  location  only 


Decreasing  K  is  good 
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The  Configuration-Dependent  Component 


MM,  v 


l  m  y=i 

Approximations  for  Z»Zr : 


-1 


•  System  I  factor  a 

1 

\6L2sL2 

- 1 - 

0 

— o-o-o-o — > 
r-1  x 

4r  I 

•  System  II  factor  5 

1 

2L4 

-LI2 

+L/2 

A/ote:  System  II  ranging  performance  depends  on  L  and 
this  component  is  independent  of  sub-aperture  length,  Ls 
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Performance  Vs.  Array  Length,  L  IJ3^, 

For  Small,  Medium,  and  Large  Ls  (with  Approximations  dashed) 


System  I.  Dual  Aperture 


System  II.  Triple  Aperture 


Form  »  Appronmabon  • 


6  Overlapping  Curves: 
Actual  and  Approximations  for 
Small,  Medium,  and  Large  Ls 

«0  SO  80  TO  80 

Length,  L  Length,  L 

•  Approximations  hold  for  L»LS 

•  System  I  depends  on  sub-aperture  length,  Ls 

•  Range  error  variance  decreases  with  increased  L  as  shown 
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Ranging  Performance  Comparison 


Range  Error  Variances  for  : 

•  System  I,  0/(i?)oc — j 


16 L2sL2 


VP 

I  Ls 


System  II)<Tj(^)«  — 


2D 


o 

-0-0*00- 


L, 


-L/2 


+L/2 


>  Take  the  square  root  of  the  ratio  of  error  variances: 


•  Relative  RMS  range  error  is  linear  with 
the  ratio  of  L  to  Ls 

•  RMS  error  expected  from  System  /  is 
35  times  that  of  System  II  for  L  =  100  Ls 

System  II  performs  better  than  System  I  for  all  cases  studied 
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°7 

f\ 

°7/ 

(*) 

V2  L_ 

~4~~n 


=  0.35— 
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Effect  of  Additional  Sensors 


Two  M-Sensor  Configurations  Relative  RMS  Range  Error 


Relative  results  are  insensitive  to  the  number  of  sensors:  M 
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Conclusions 


MAX/SEA 


***HO*T 

•  For  the  two  sonar  systems  considered  with  the 
same  aperture  (length:  L)  and  the  same  number  of 
sensors:  M,  the  relative  ranging  performance  was 
calculated.  The  theory  has  various  applications. 

•  For  examples  considered,  improvement  factors  of 
35  to  40  in  RMS  range  error  were  observed. 

•  Any  perceived  cost  advantage  gained  by  installing 
sensors  at  only  two  sites  instead  of  three  must  be 
weighed  against  the  significantly  improved  ranging 
performance  of  a  three-site  system  unless  the  two- 
site  system  meets  the  required  performance. 
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